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Abstract. Ball throwing on Euclidean spaces has been considered for a while. A 
suitable renormalization leads to a fractional Brownian motion as limit object. In this 
paper we investigate ball throwing on spheres. A different behavior is exhibited: we 
still get a Gaussian limit but which is no longer a fractional Brownian motion. However 
the limit is locally self-similar when the self-similarity index H is less than 1/2. 
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1. Introduction 

Random balls models have been studied for a long time and are known as germ- 
grain models, shot-noise, or micropulses. The common point of those models consists 
in throwing balls that eventually overlap at random in an n-dimensional space. Many 
random phenomena can be modelized through this procedure and many application fields 
are concerned: Internet traffic in dimension one, communication network or imaging in 
dimension two, biology or materials sciences in dimension three. A pioneer work is due 
to Wicksell [21] with the study of corpuscles. The literature on germ- grain models deals 
with two main axes. Either the research focuses on the geometrical or morphological 
aspect of the union of random balls (see [20] or [2T] and references therein), or it is 
interested in the number of balls covering each point. This second point of view is 
currently known as shot-noise or spot-noise (see [6] for existence). In dimension three, 
the shot-noise process is a natural candidate for modeling porous or granular media, and 
more generally heterogeneous media with irregularities at any scale. The idea is to build 
a microscopic model which yields a macroscopic field with self-similar properties. The 
same idea is expected in dimension one for Internet traffic for instance [25J. A usual way 
for catching self-similarity is to deal with scaling limits. Roughly speaking, the balls are 
dilated with a scaling parameter A and one lets A go either to or to infinity. We quote 
for instance [1] and [T2] for the case A — ► + , [TT] and [2J for the case A — > +oo, [3] and 
[3] where both cases are considered. 

In the present paper, we follow a procedure which is similar to |2J and [3J. Let us 
describe it precisely. A collection of random balls in M. n whose centers and radii are 
chosen according to a random Poisson measure on R n x R + is considered. The Poisson 
intensity is prescribed as follows 

-n-l+2H 
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for some real parameter H. Since the Lebesgue measure dx is invariant with respect 
to isometry, so is the random balls model, and so will be any (eventual) limit. As the 
distribution of the radii follows a homogeneous distribution, a self-similar scaling limit 
may be expected. Indeed, with additional technical conditions, the scaling limits of 
such random balls models are isometry invariant self-similar Gaussian fields. The self- 
similarity index depends on the parameter H. When < H < 1/2, the Gaussian field 
is nothing but the well-known fractional Brownian motion [16} [18| [19] . 

Manifold indexed fields that share properties with Euclidean self-similar fields have 
been and still are extensively studied (e.g. [HI El UHl H3 [T71 [22l E3]). In this paper, 
we wonder what happens when balls are thrown onto a sphere, and no longer onto a 
Euclidan space. More precisely, is there a scaling limit of random balls models, and, 
when it exists, is this scaling limit a fractional Brownian field indexed by the surface for 
< H < 1/2? 

The random field is still obtained by throwing overlapping balls in a Poissonian way. 
The Poisson intensity is chosen as follows 

^(dx,dr) = f(r)a(dx)dr. 

The Lebesgue measure dx has been replaced by the surface measure a(dx). The function 
/, that manages the distribution of the radii, is still equivalent to r - n ~ 1 + 2H ; a t least 
for small r, where n stands for the surface dimension. It turns out that the results are 
completely different in the two cases (Euclidean, spherical). In the spherical case, there 
is a Gaussian scaling limit for any H. But it is no longer a fractional Brownian field, as 
defined by [p3]. We then investigate the local behavior, in the tangent bundle, of this 
scaling limit in the spirit of local self-similarity [H [71 [T5] . It is locally asymptotically 
self-similar, with a Euclidean fractional Brownian field as tangent field. Our microscopic 
model has led to a local self-similar macroscopic model. 

The paper is organized as follows. In Section 2, the spherical model is introduced and 
we prove the existence of a scaling limit. In Section 3, we study the locally self-similar 
property of the asymptotic field. Section 4 is devoted to a comparative analysis between 
Euclidean and spherical cases. Eventually, some technical computations are presented 
in the Appendix. 

2. Scaling limit 
We work on S n the n-dimensional unit sphere, n > 1: 

§n = {(Xi)l<i< n+ 1 € 5 Yl X i =1 }- 

l<j<n+l 

2.1. Spherical caps. For x,y G S n , let d(x,y) denote the distance between x and y on 
S n , i.e. the non-oriented angle between Ox and Oy where O denotes the origin of R n+1 . 
For r > 0, B(x,r) denotes the closed ball on S centered at x with radius r: 

B(x,r) = {y G §„ ; d(x,y) < r}. 

Let us notice that for r < tt, B(x,r) is a spherical cap on the unit sphere § n , centered 
at x with opening angle r and that for r > ir, B(x, r) = S n . 
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Denoting cr(dx) the surface measure on S n , we prescribe 4>(r) as the surface of any ball 
on S n with radius r, 

4>{r) := a(B(x,r)) , x £ S n ,r > . 

We also introduce the following function defined for z and z', two points in S n and 

r £ R+, 

(1) V(z,z',r):= ld(x,z)<r ld(x,z')<r °"( d2 • 

Actually ^(z, z', r) denotes the surface measure of the set of all points in S n that belong 
to both balls B(z, r) and B(z', r). Clearly z', r) only depends on the distance d(z, z') 
between z and z' . We write 

(2) i>(d(z,z'),r) = t!(z,z',r) 

and note that it satisfies the following: V(u, r) £ [0, n] x R + 

• < ip(u,r) < a(§ n ) A<j>(r) 

• if r < u/2 then ip(u, r) = and if r > ir then ^("U) r ) = 0"(§n) 

• ^(0, r) = 0(r) ~ cr n as r — > + . 

In what follows we consider a family of balls B(Xj,Rj) generated at random, following 
a strategy described in the next section. 

2.2. Poisson point process. We consider a Poisson point process (Xj, Rj)j in S„xl + , 
or equivalently iV(dx,dr) a Poisson random measure on S n x R + , with intensity 

j/(dx, dr) = /(r) <r(dx) dr 

where / satisfies the following assumption A(H) for some H > 0: 

• supp(f) C [0, vr) 

• / is bounded on any compact subset of (0, it) 

• /(r) ~ r~ n - 1+2H as r 0+. 

Remarks : 

1) The first condition ensures that no balls of radius Rj on the sphere self- intersect. 

2) Since (f>(r) ~ cr n , r — > + , the last condition implies that J R+ <p(r)f(r)dr < +oo, 
which means that the mean surface, with respect to /, of the balls B(Xj, Rj) is finite. 

2.3. Random field. Let Ai denote the space of signed measures [i on § n with finite 
total variation |//|(S n ), with the total variation measure of fi. For any fi £ Ai, we 
define 

(3) X(p)= J n(B(x,r)) N(dx,dr) . 
Note that the stochastic integral in ([3]) is well defined since 

/ \n(B(x,r))\ f(r)a(dx)dr < f f f l d{x , y)<r f(r) a(dx) |yu|(dy) dr 

JS n xM.+ JS n Js n JR+ 

= \fi\(S n ) U<f>(r)f(r)dr\ < +00 . 



4 ANNE ESTRADE AND JACQUES ISTAS 

In the particular case where \i is a Dirac measure S z for some point z G S> n we simply 
denote 

(4) X(z) = X(5 Z ) = f l B (x,r)(z) N(dx,dr) . 

i§ n xM+ 

The pointwise field {X(z); z G §> n } corresponds to the number of random balls (Xj, Rj) 
covering each point of § n . Each random variable X(z) has a Poisson distribution with 
mean J" R+ 4>{r)f(r) dr. 
Furthermore for any G «M, 

E(X(/z)) = M§„) ( / </>(r)f(r)dr 

V./R+ 

and 

Var(X(/x)) = / n(B(x,r)ff(r)a{dx)dr G (0, +oo] . 

2.4. Key lemma. For H > 0, we would like to compute the integral 

f n(B(x, r)) 2 r - n ~ 1+2H a(dx) dr 
JS n xR+ 

which is a candidate for the variance of an eventually scaling limit. We first introduce 
A4 H the set of measures for which the above integral does converge: 

M H = M if 2H < n ; M H = {n G M ; /i(S„) = 0} if 2if > n . 

The following lemma deals with the function tp prescribed by ([2]). 

Lemma 2.1. Let H > u>i£/i 2i7 7^ n. VFe introduce 

ip( H ) =ipifO<2H<n; ip {H) = ip - a{§ n ) if 2H > n . 

Then for all u G [0, 7r], 

/ \^ H \u,r)\r- n ~ 1+2H dr < +00. 

Furthermore, denoting 

(5) = / V (H) (-"^)^" n ~ 1+2H dr 
/or any u in [0,7r], we /iaue for all [i G .M^ ; 

0< f fi(B(x,r)) 2 r- n - 1+2H a(dx)dr = [ K H (d(z, z')) fi(dz)fi(dz') < +00. 



Remark 2.2. 

i,) For x,y in S n , the difference of Dirac measures 5 X — 5 y belongs to M. H for any H. 

2) In the case 2H > n, since any fi G Ai H is centered, the rhs integral is not changed 
when a constant is added to the kernel Kh- 

3) This lemma proves that the kernel Kh defines a covariance function on A4 H . 
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Proof. Using the properties of ip, we get in the case < 2H < n, 
0< f 4>(u,r)r- n ~ 1+2H dr < [ 4>(r) r - n - 1+2H dr + a(S n ) [ r - n - 1+2H dr 

< + OO . 

In the same vein, in the case 2H > n, we get 

0< / K§ n ) - 1>(u,r)) r~ n - 1+2H dr < a(S n )[ r ~n-i+2H dr 

< +OO . 

We have just established that there exists a finite constant Ch such that 

(6) Vu€[0,tt], f \^ H) (u,r)\r- n - 1+2H dr<C H - 

The first statement is proved. 

Let us denote for /j, G M H 

I H (») = f fi(B{x,r)) 2 r- n - 1+2H a{dx)dr 

and start with proving that Ih{h) is finite. We will essentially use Fubini's Theorem in 
the following lines. 

W =1(1 KB(x,r)) 2 a(dx))r~ n - 1+2H dr 



( [ *(z, z', r) fi(dz)fi(dz')) r -«"i+2^ dr . 

Since fi € A4 H is centered in the case 2H > n, one can change i/j into ip( H ) within the 
previous integral. Then 

Ih(jj) < f ([ |^ (H) (^,^) ) r)||Ml(d^)|/i|(d^)V _n " 1+2Hdr 

< f ([ \^ H \d(z,z'),r)\r- n ~ 1+2H dr\\^\{dzM(dz ! ) 

< C H \fi\ (§n) 2 < +00 . 



Following the same lines (except for the last one) without the "| |" allows the computation 
of Ih(i^) and concludes the proof. □ 

An explicit value for the kernel Kjj is available starting from its definition. The point 
is to compute ip^ H '. We give in the Appendix a recurrence formula for ip( H \ based on 
the dimension n of the unit sphere § n (see Lemma l4.ip . 
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2.5. Scaling. Let p > and A be any positive function on (0,+oo). We consider the 
scaled Poisson measure N p obtained from the original Poisson measure N by taking the 
image under the map (x,r) G S x M + (x,pr) and multiplying the intensity by X(p). 
Hence N p is still a Poisson random measure with intensity 

v p (dx, dr) = \(p)p~ l f(p~ l r) a(dx) dr . 

We also introduce the scaled random field X p defined on M. by 

(7) X P (j*)= I p(B(x,r)) N p (dx,dr) . 

JS n xR+ 

Theorem 2.3. Let H > with 2H ^ n and let f satisfy A(H). For all positive 
functions A such that \{p)p n ~ 2H — > +oo, the limit 

( Xp (p)-E(x p (p)) m h\ m {Wh ^,^M H } 



holds in the sense of finite dimensional distributions of the random functionals. Here 
Wh is the centered Gaussian random linear functional on J\4 H with 

(8) Cov(W H (p),W H (v))= j K H (d(z,z'))p(dz)v(dz'), 



where Kh is the kernel introduced in Lemma \2.1\ 

The theorem can be rephrased in term of the pointwise field {X(z); z G S n } defined 
in (jH). 

Corollary 2.4. Let H > with 2H ^ n and let f satisfy A(H). For all positive 
functions A such that X(p)p n ~ 2H — > +oo, 

• if < 2H < n then 

{^^■■^)^^ 

where Wh is the centered Gaussian random field on § n with 

Cov(W H (z),W H (z')) = K H {d(z,z')) . 

• if 2H > n then for any fixed point zq G S n , 

f X p (z) - X p (z ) \ fdd 

where Wh, zo ^ s the centered Gaussian random field on § n with 

Cov(W HiZ0 (z),Wh, Z0 (z')) = K H (d(z,z')) - K H (d(z,z )) - K H (d(z',z )) + K H (0) . 
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Proof, of Theorem 

Let us denote n(p) := \J \(p)p n ~ 2H . The characteristic function of the normalized field 
(X p (.) - E(X p (.))) /n(p) is then given by 

(9) E ( mp -wm y = „ (j^ arM ^ 



where 

(10) G p (x, r) = [e l ^T- - 1 - ^"^'" j Kp)p~ 1 f{P~ l r) 

We will make use of Lebesgue's Theorem in order to get the limit of L xR+ G p (x, r) dra(dx) 
as p — > +00. 

/ ■/J 1 (B(x,r)) (~D( \\ \ 

On the one hand, n{p) tends to +00 so that ( e n < p > — 1 — i — w (pj j behaves like 



( n(B(x,r)) x 2 
n(p) 



. Together with the assumption A(H), it yields the following asymptotic. 



For all (x,r) e§ n xl+, 

(11) G p (x,r) — > - l -p{B{x,r)?r- n - l + 2H . 

On the other hand, since ^^(p)^ — \p\(^n) f° r P l ar § e enough, we note that there 
exists some positive constant K such that for all x,r,p, 



i ^B(x, r )) a(B(x,r)) 



<R f p(B(x,r)) 

V n (p) 



n(p) 

Therefore 

|G p (;r,r)| < Kp(B(x,r)) 2 p- n - 1+2H f(p- 1 r). 

There exists C > such that for all rEl + ,/(r)< Cr- n ~ l+2H . Then we get 

(12) \G p (x,r)\ < KC p(B(x,r)) 2 r- n - 1+2H 

where the right hand side is integrable on § n x M + by Lemma 12.11 
Applying Lebesgue's Theorem yields 

1 

/s„xi+ p^+oo 2 j §nX 

1 

~2 

Hence (X p (p) — ~E(X p (p)))/n(p) converges in distribution to the centered Gaussian 
random variable W(p) whose variance is equal to 

E(W(p) 2 )=C I K H (d(z,z'))p(dz)p(dz') . 
By linearity, the covariance of W satisfies flHJ) . 

□ 



G p (x,r)a{dx)dr — ► -- / p(B(x, r)) 2 r - n - 1+2H a(dx) dr 

(>->+°° 2 7s n xR+ 

K H {d{z,z'))p{dz)p(dz') . 
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Remark 2.5. 

1) The pointwise limit field {Wh(z);z G §> n } in Corollary \2.J\ is stationary, i.e. its 
distribution is invariant under the isometry group of S n , whereas the increments of 
{Wh,z ( z )', z £ §n} ar & distribution invariant under the group of all isometries of S n 
which keep the point zq invariant. 

2) When < H < 1/2 the Gaussian field Wh does not coincide with the field introduced 
in [13] as the spherical fractional Brownian motion on § n . 

Indeed, let us have a look at the case n = 1, it is easy to obtain the following piecewise 
expression for ip = ij}^: \/(u,r) € [0, tt] x R + , 

ipi(u,r) = for < r < u/2 

= 2r — u for u/2 < r < ir — u/2 
= 4r — 2ir for ir — u/2 < r < ir 
= 2ir for r > 7r 

and to compute 

Kh{u) = gft _ \ H)22H (2(2H)™ - - (2n - uf H ) , 

Actually, we compute the variance of the increments of Wh 
E(W h (z)-W h (z')) 2 = 2K H (0)-2K H (d(z,z')) 

= H{l-2H)2^ [d2H{z ' Z ' ] + {27T ~ Z ' )?H ~ {27T)2H] - 

The spherical fractional Brownain motion Bh, introduced in [13], satisfies 

E(B H (z) - B H {z')f = d 2H (z,z'). 

Even up to a constant, processes Wh and Bh are clearly different. The Euclidean situa- 
tion is therefore different. Indeed, [3], the variance of the increments of the corresponding 
field Wh is proportional to \z — z'\ 2H . 

3. Local self-similar behavior 

We wonder whether the limit field Wh obtained in the previous section satisfies a 
local asymptotic self-similar (lass) property. More precisely we will let a "dilation" of 
order e act on Wh near a fixed point A in S n and as in [15], up to a renormalization 
factor, we look for an asymptotic behavior as e goes to 0. An //-self-similar tangent 
field Th is expected. Recall that Wh is defined on a subspace A4 H of measures on § n , 
so that Th will be defined on a subspace of measures on the tangent space T^n of S n . 

3.1. Dilation. Let us fix a point A in § n and consider T^Sn the tangent space of § n at 
A. It can be identified with M. n and A with the null vector of W 1 . 

Let 1 < 5 < tt. The exponential map at point A, denoted by exp, is a diffeomorphism 

o 

between the Euclidean ball {y £ W 1 , \\y\\ < 5} and B(A,5) C S n , where ||.|| denotes the 
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o 

Euclidean norm in M n and B(A, 5) the open ball with center A and radius 5 in S n . 
Furthermore for all y,y' £ M n such that ||y||, < 5, 

d(A,expy) = \\y\\ and d(exp y, exp y) < \\y-y'\\ . 

We refer to [TO] for precisions on the exponential map. 

Let r be a signed measure on W 1 . We define the dilated measure t £ by 

VBGB(M n ) r B (B) = t(5/e) 

O 

and then map it by the application exp, defining the measure \x e = exp* r e on B(A, 5) 
by 

(13) VC£B(B(A,5)) fi £ (C) = exp* t £ (C) = r^exp-^C)). 

We then consider the measure /x E as a measure on the whole sphere S n with support 

o 

included in B(A, 6). 

At last, we define the dilation of Wh within a "neighborhood of A" by the following pro- 
cedure. For any finite measure r on R n , we consider fx £ = exp* t £ as defined by (|13p and 
compute Wjj(fie)- We will establish the convergence in distribution of e~ H Wn(exp* t £ ) 
for any r in an appropriate space of measures on W 1 . Since Wh{^e) is Gaussian, we will 
focus on its variance. 

3.2. Asymptotic of the kernel Kh- For < H < 1/2, we already mentioned that 
the kernel Kh(0) — Kh{u) is not proportional to u 2H . As a consequence, one cannot 
expect Wh to be self-similar. Nevertheless, as we are looking for an asymptotic local 
self-similarity, only the behavior of Kh near zero is relevant. Actually we will establish 
that, up to a constant, Kh(0) — Kh{u) behaves like u 2H when u — > + . 

Lemma 3.1. Let < H < 1/2. The kernel Kh defined by (jSj) satisfies 

K H (u) =Ki- K 2 u 2H + o(u 2H ), u^0 + 

where K\ = Kh(0) and K 2 are nonnegative constants. 

Proof. Let us state that the assumption H < 1/2 implies H < n/2 so that in that case 
Kh is prescribed by 

K H {u) = [ ij(u,r)r~ n - l+2H dr , u G [0,tt] . 
Jr+ 

We note that Kh(0) < +oo since ^(0, r) ~ cr n as r — > + and ^(0, r) = cr(S n ) for r > n. 
Then, subtracting X^(0) and remarking that ip(0,r) = ip(u,r) = <r(S n ) for r > ir, we 
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write 

K H (0)-K H (u) = / (V>(0,r) - ifj(u,r))r- n - 1+2H dr 
Jo 

6 



(V>(0, r) - i/)(u, r)) r - n - 1+2H dr 

f'TT 

+ (t/j(0,r) -ip(u,r))r~ n - 1+2H dr 



where we recall that 5 £ (l,7r) is such that the exponential map is a diffeormorphism 
between {||y|| < 6} C R n and B(A,S) C §„. 

The second term is of order u, and therefore is negligible with respect to u 2H , since tp is 
clearly Lipschitz on the compact interval [5, ir] . 

We now focus on the first term. Performing the change of variable r i— > r/u, we write it 
as 

rS r 
(^(0, r) - r)) r - n - 1+2H dr = u 2H / A(u, r) r - n - 1+2H dr , 



where 

A(u,r) := l ur<s u~ n (ip(u,ur) -ip(0,ur)) . 

Their only remains to prove that f R+ A(u,r) r~ n ~ 1+2H dr admits a finite limit K.% as 
u — ► + . We will use Lebesgue's Theorem and start with establishing the simple con- 
vergence of A(u,r) for any given r G M + . 

We fix a unit vector v in R n and a point A' = expv in S n . We then consider for 
any u £ (0, 5), the point := exp(uv) G § n located on the geodesic between A and A' 
such that ^44) = [|itv|| = u. We can then use (UJ and ([2]) to write 

V>0,.) = $(A,A' U ,.) = / l d (A, z )<. ld(A', z )<. da(z) 



if>(0, .) = A, .) = / l d(Az)< . d(r(z) 
in order to express A(u, r) as 

A(n,r) = l ur<5 U~ n / l rf (A, 2 )< U r l<f(^,*)>ur dcr(» • 

o o 

Since ur < 5 the above integral runs on B(A,ur) C 5) and we can perform the 

exponential change of variable to get 

A(u,r) = l ur< g U~ n / l|| y ||<ur ld(exp(uv),exp(y))>urd0'( ex P(j/)) 
^-ur<8 I ^-\\y\\<r ^-d(exp(wv),exp(uj/))>ur &\ u y)dy . 

In the last integral, the image by exp of the surface measure d<r(exp(y)) is written as 
&{y)dy where dy denotes the Lebesgue measure on W 1 . 
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We use the fact that d(exp(ux), exp(ux')) ~ u\\x — x'\\ as u — > + to get the following 
limit for the integrand 

l<i(exp(uv),exp(Mj/))<ur ^( u lj) * ^||v— y||>r • 

Since the integrand is clearly dominated by 

IHloo := sup{5-(y) , ||y|| < 5} , 
Lebesgue's Theorem yields for all r £ R + , 

A(u,r) — > cr(0) / l[i y [| <r l|i v _ y [| >r dy . 

We recall that d(exp x,exp x') < \\x — x'\\ for all x,x' G with norm less than 5. 
Therefore for all u, 

A(u,r) < \\a\\oo / l\\ y \\ <r l\\ v - y \\ >r dy 
where the right hand side belongs to L 1 (IR + , r _n_1+2 ^ dr) (see [2] Lemma A. 2). 
Using Lebesgue's Theorem for the last time, we obtain 



A(u,r)r- n - 1+2H dr — ► K 2 



where 



(0) / l N | <r l|| v _,,|| >r r-' t - 1+2 ^dydr G (0,+oo) . 

JR"xl+ 



□ 



Let us remark that the proof makes it clear that the case H > 1/2 is dramatically 
different. The kernel Kh(0) — K — H(u) behaves like u near zero and looses its 2H 
power. 



3.3. Main result. 

Let < H < 1/2. We consider the following space of measures on T^n — M n 
9Jl H = {measures r on M. n with finite total variation such that 



and 



/ \\x-x'\\ 2H \T\(dx)\T\(dx') <+00} . 

JR n xR n 



For any measure r G dJt H , we compute the variance of Wh(jj, e ) where /i E = exp* r £ is 
defined by (TBI) . 

By Lemma |2.1( since fj, £ belongs to Ai = M H in the case H < 1/2, 

var^O^)) = f K H (d(z,z'))fi £ (dz)fi £ (dz'). 
JB(A,5)xB(A,S) 
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Performing an exponential change of variable followed by a dilation in M n , we get 
var(Wff(/i e )) = / 1 

||y||<<5 l|[j/'||<5 

«"xM" 

HxWkS/s Hx'WkS/s K H (d(exp(£x), exp(ex')))r(dx)r(dx / ). 



Denoting Kh(u) = Kh(u) — Kh(0), 

vaT(W H (n £ )) = / l\\ x \\ <s/E l\\ x/ \\ <s/£ K H (d(exp(ex),exp(ex')))T(dx)T(dx') 

JR n xR n 

+ K H (0) t({\\x\\ <5/e}) 2 . 
Let us admit for a while that 

r({ M <*/ e} )» _ ^ 

Then, applying Lebesgue's Theorem with the convergence argument on Kh obtained in 
Lemma yields 

(15) Var( ^ )} — -kJ \\x-x'f H r{dx)r(dx'). 

Let us now establish (|14p where we recall that r is any measure in 9Jl H . In particular, 
the total mass of r is zero so that 

r({\\x\\ < 5/e}) r({\\x\\>S/e}) 



J 



£ H l\\ x \\ >S / £ T(dx) 



For any fixed x £ M n , e ^ l||x||><5/£ * s zero wrien e is small enough. Moreover e H l\\ x \\>s/e 
is dominated by 5~ H \\x\\ which belongs to L 1 (M™, |r|(dx)) since r belongs to 9Jt H . 
Lebesgue's Theorem applies once more. 



We deduce from asymptotic (|15p the following theorem. 
Theorem 3.2. Let < H < 1/2. The limit 

W#(exp*T e ) fdd 



/zoZefe for all r G 971^, in the sense of finite dimensional distributions of the random 
functionals. Here Tj{ is the centered Gaussian random linear functional on 9Jt H with 

(16) Cov(T h {t),T h (t')) =-K 2 I ||x-x , || 2 ^r(dx)r , (dx / ), 

JM. n xR n 

As for Theorem 1 2. 3 \ Theorem 13.21 can be rephrased in terms of pointwise fields. Indeed, 
&x — So belongs to 9Jl H for all x in W 1 . Let us apply Theorem 13.21 with r = 5 X — do- 
Then Th(5 x — 5o) has the covariance 

Cov (T H (5 X - 5o),T H (5 x , - So)) = K 2 (\\x\\ 2H + \\x'\\ 2H - \\x - x'\\ 2H ) , 
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and the field {Th(S x — <5o) ', x £ R n } is a Euclidean fractional Brownian field. 



4. Comparative analysis 



In this section, we aim to discuss the differences and the analogies between the Eu- 
clidean and the spherical case. 

Let us first be concerned with the existence of a scaling limit random field. The variance 
of this limit field should be 



where M n is the n-dimensional corresponding surface with its surface measure a. When 
speaking of the Euclidean case M n = W 1 we refer to [3]. In the present paper, we 
studied the case MI n — S n . Moreover, in this discussion, the hyperbolic case M n — EI n — 
{(xi)i<i< n+1 £ R n+1 ; x 2 n+l - J2i<i<n x i = 1 > x n+i > 1} is evoked. 

In the Euclidean case, the random fields are defined on the space of measures with 
vanishing total mass. So let us first consider measures fi such that /x(M n ) = 0. Hence, 
whatever the surface M n , the integral V involves the integral of the surface of the 
symmetric difference between two balls of same radius r. As r goes to infinity, three 
different behaviors emerge. 

• M n = S n : this surface vanishes 

• M n = W 1 : the order of magnitude of this surface is r n ~ l 

• M n = M n : the surface grows exponentially 

The consequences are the following. 

• M n = S n : any positive H is admissible 

• M n = R n : the range of admissible H is (0, 1/2) 

• M n = H n : no H is admissible 

In the Euclidean case, the restriction fj,(W l ) = is mandatory whereas it is unnecessary 
in the spherical case for H < n/2. Indeed the integral V is clearly convergent. 

Let us now discuss the (local) self-similarity of the limit field. Of course, we no longer 
consider the hyberbolic case. 

• M ra = M. n : dilating a ball is a homogeneous operation. Therefore, the limit field 
is self-similar. 

• M n = S n : dilation is no longer homogeneous. Only local self-similarity can 
be expected. The natural framework of this local self-similarity is the tangent 
bundle, where the situation is Euclidean. Hence we have to come back to the 
restricting condition H < 1/2. 



Recall that the functions ip n 's are defined by ([I]) and (|2|) 

*p n (u,r) = ty n (M,M',r) = / l d{M ,N)<r U(M',N)< r da n (N) , (u,r) G [0, vr] x R 



V 




Appendix 



Recurrence formula for the ^» n 's. 
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for any pair (M, M') in S n such that d(M,M') = u. Here a n stands for the surface 
measure on S n . 

Lemma 4.1. The family of functions tp n ,n > 2 satisfies the following recursion: 
V(«, r) G [0, tt] x R+ , 



t (u,r) = y (1 - a 2 ) w/2 V>»-i ^,arccos (^ JpL =^ da . 



Proof. An arbitrary point of S n is parameterized either in Cartesian coordinates, (xj)i<j< n+ i, 
or in spherical ones 

with 







l<i<n 


€ [O,^)™- 1 x [0,2vr) 


Xl 




COS 


01 


X2 




sin i 


f>l COS 02 


x 3 




sini 


sin cf>2 cos 03 



x n = sin 0i sin 02 • • • sin0 n _i cos n 
x n+ i = sin 4>i sin <p 2 ■ ■ ■ sin n _ i sin n 

Let M be the point (0i)i<j< n = (0, . . . , 0). One can write the ball B n (M, r) of radius 
r, which is a spherical cap on S n with opening angle r as follows, 

B n (M,r) = {(&)i<i<„ G S n ; 01 < r} 

or in Cartesian coordinates 

B n (M,r) = {{xi)i<i< n+l £ S n ; xi > cosr} . 

Let a 6 (—1,1) and let P a be the hyperplane of W l+1 defined by x n+ \ = a. Let us 
consider the intersection P a n B n (M, r). 



• If 1 - a 2 < cos 2 r then P a n B n (M, r) 



• If 1 — a 2 > cos 2 r then 
P a n B n (M,r) = {(xj)i<j< n+ i € S n ; xi > cosr and x n+ i = a} 

= {{xi)\<i< n G M n ; xi > cosr and ^ xf = 1 — a 2 } x {a}. 

l<j<n 

In other words, denoting S n _i(i2) the (n — l)-dimensional sphere of radius R, 
P a n B n (M,r) = B n _ ljVT ^(M(a),r(a)) x {a} 



where B n l ^^p(M(a),r(a)) is the spherical cap on S n _i(\/l — a 2 ), centered 
at M(a) = (Vl — a 2 , 0, . . . , 0) and with opening angle r(a) = arccos 



Let now M' be defined in spherical coordinates by (0j)i<j< n = (u, 0, ...,0), so that 
d(M, M') = u. The intersection P a n B n (M',r) is the map of P a n B n (M,r) by the 
rotation of angle u and center C in the plane x$ = . . . = x n+ \ = 0. So 
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• if 1 - a 2 < cos 2 r then P a n B n (M' , r) = 0. 

• if 1 — Xq > cos 2 r then 

P a nB n (M',r) = B n _ ltVT ^(M'{a),r(a)) x {a} 

where the (n — l)-dimensional spherical cap B nl y / 1 _ a 2 (M' (a),r(a)) is now cen- 
tered at M'(a) = (vl — a 2 cos u, \/i — a 2 sinu, 0, . . . , 0). 

We define tp n -i,R(u, r) as the intersection surface of two spherical caps on S n _i(i2), 
whose centers are at a distance Ru and with the same opening angle r. 
By homogeneity, this leads to 

ip n - hR {u,r) = RJ 1 ^ 1 ip n _ lt i(u,r) = R n ip n -i{u, r) . 

The surface measure a n of §„ can be written as 



a 2 do 



da n (xi, ...,x n ,a) = \J\-a 2 da n ^ ly/T ^j(xi, ...,x n ) x da 

where cT n -i,_R is the surface measure of S n _i(i?). 
We then obtain 

/ ( COS T \ I — 
^ ll-a2>co S 2r V' rl _i iV T=^( n ' arCCOS \ —B _ fl2 )) V 1 

sinr / cosr \ 

(VI - a 2 )™ Vn-i(«,arccos )da, 

-sinr VVl — a / 

and Lemma |4. II is proved. □ 
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